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SUMMARY 


[ This  report  describes  certain  results  obtained  in  the  course  of 

i a study  of  methods  for  the  analysis  of  the  high  frequency  vibration  of  struc- 

i tural  systems.  The  main  problems  which  arise  are  the  need  for  a fine  enough 

(• 

[ spatial  discretization  to  transmit  accurately  the  high  frequency  vibrations, 

; the  need  for  a suitable  time  integration  algorithm  that  is  accurate  for  the 

high  frequency  range  and,  finally,  the  need  for  an  accurate  representation  of 
r the  damping  characteristics  of  the  structure. 

To  gain  insight  into  these  problem  areas,  a series  of  idealized 
models  are  considered  which  incorporate  the  characteristics  of  an  embedment 
media,  a structure  and  internal  equipment,  the  equipment  having  natural  fre- 
quencies which  are  high  compared  with  the  fundamental  frequency  of  the  struc- 
ture. The  purpose  of  the  analysis  is  to  compare  exact  solutions  to  the 
; idealized  model  with  appropriate  solutions,  which  may  be  obtained  by  standard 

structural  dynamics  techniques  to  give  an  understanding  of  the  degree  of  re- 
finement needed  in  the  modeling  and  the  numerical  technique.  Certain  results 
on  the  vibration  of  equipment  have  been  obtained  which  appear  to  be  new. 

These  techniques  were  applied  to  a structure-equipment  system 
using  more  realistic  models  that  are  representative  of  silo-like  protective 
structures  and  it  was  shown  that  essentially  the  same  type  of  results  developed. 
This  shows  that  the  phenomenon  is  not  confined  to  the  elementary  idealized 
model.  Excitations  simulating,  ground  shock,  air  shock  or  shaker  loading  were 
examined. 

An  alternative  simple  modeling  of  structural  damping  is  intro- 
duced which  has  much  the  same  degree  of  simplicity  as  the  usual  damping  model 
but  leads  to  improved  estimates  of  structural  response  at  high  frequencies. 

The  model  is  based  on  the  well  known  Maxwell  model  of  viscoelasticity  and  in 
a one  degree-of-freedom  system  requires  only  one  damping  parameter. 

For  sinusoidal  loading  of  a single  degree-of-freedom  system  the 
response  is  that  of  a conventional  damping  around  resonance  but  differs  at 


high  frequencies.  In  a multi  degree-of-freedom  system  it  is  possible  to  use 
modal  damping  facto.s.  It  is  shown  tha*"  algorithms  for  the  direct  numerical 
integration  of  the  matrix  equations  of  motion  of  the  multi  degree-of-freedom 
system  can  be  easily  adapted  to  this  type  of  modeling. 
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SECTION  1 


INTRODUCTION 

Conventional  methods  for  the  dynamic  analysis  of  large  structural 
systems  are  designed  to  be  accurate  in  the  computation  of  the  response  of  the 
system  for  the  lowest  modes.  In  most  applications,  as  for  example  In  earth- 
quake engineering  design,  only  the  low-mode  response  is  of  Interest  and  in 
these  cases  the  use  of  implicit  unconditionally  stable  integration  algorithms 
is  preferred.  In  addition  to  being  unconditionally  stable,  when  low-mode 
accuracy  is  needed  it  is  advantageous  for  an  algorithm  to  possess  some  form 
of  numerical  dissipation  to  damp-out  any  spurious  participation  from  higher 
modes.  Commonly  used  algorithms  in  structural  dynamics  which  possess  these 
characteristics  are  Houbolt's  method  (Reference  1),  Wilson's  6 method  (Refer- 
ence 2)  and  the  Newmark  family  of  methods  (Reference  3).  In  earthquake 
engineering  these  methods  are  justified  on  the  basis  that  the  higher-mode 
response  does  not  affect  the  structural  design.  Another  reason  for  suppress- 
ing the  higher  frequencies  is  that  size  limitation  forces  a discretization 
which  does  not  reproduce  the  higher  modes  accurately. 

There  are  many  situations  in  which  the  high  frequency  response  of 
a structure  might  be  of  interest,  for  example  it  may  be  necessary  to  estimate 
the  accelerations  induced  in  non-structural  elements  in  buildings,  such  as 
equipment  or  ordinance  and  it  may  also  be  necessary  to  design  mountings  for 
sensitive  equipment.  Most  computer  programs  used  by  structural  designers 
utilize  a fairly  crude  modeling  of  the  system  considering,  for  example,  only 
three  degrees  of  freedom  per  story  and  modeling  damping  using  either  viscous 
or  hysteretic  models.  These  techniques  are  accurate  in  a typical  structural 
frame  up  to  frequencies  around  10  Hz.  Equipment  design  and  design  of  non- 
structural  elements  may  require  techniques  which  are  accurate  to  around  200  Hz. 
The  conventional  model  will  thus  operate  as  an  analytical  filter  and  prevent 
energy  at  high  frequencies  from  being  included  in  the  motion  of  the  relevant 
non-structural  components  during  the  analysis. 

The  research  work  in  this  Report  describes  results  which  are  im- 
portant in  the  development  of  suitable  computational  techniques  to  allow  the 
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accurate  computation  at  a reasonable  cost  of  the  response  of  non-structural 
components  and  equipment  in  structures  subject  to  dynamic  loading  produced 
by  ground  and  air  shock. 


1.1.  BACKGROUND— DIFFICULTIES  WITH  CURRENT  METHODS  OF  ANALYSIS 

1.1.1.  PROBLEMS  OF  MESH  SIZE  UNIFORMITY 

The  finite  element  method  is  the  most  popular  technique  now  in 
use  for  the  dynamic  analysis  of  structures.  In  conjunction  with  large  com- 
puters the  method  offers  great  flexibility  and  accuracy.  However,  if  three- 
dimensional  structural  behavior  is  to  be  included,  an  analysis  using  this 
method  will  require  much  comptiter  time  and  considerable  expense.  To  reduce 
the  cost  of  computer  time,  many  current  programs  employ  assumptions  of  floor 
diaphragm  rigidity  and  also  neglect  inertia  at  column-floor  intersections 
which,  with  the  use  of  static  condensation,  greatly  reduce  the  number  of 
dynamic  degrees  of  freedom  to  be  considered.  These  assumptions  are  very 
useful  in  reducing  the  computer  time  of  the  analysis  but  then  also  reduce  the 
frequency  transmissibllity  of  the  model. 

Modeling  the  high  frequency  response  of  a structure  using  finite 
element  methods  requires  elements  to  be  small  and  leads  to  large  numbers  of 
elements  and  equations  of  motion  and  to  related  problems  of  storage  and 
computer  time.  Experience  in  the  use  of  codes  designed  to  solve  wave  propaga- 
tion problems  has  shown  (References  4 and  5)  that  it  is  essential  to  select  a 
mesh  which  is  as  spatially  homogeneous  as  possible  to  avoid  spurious  back 
reflections  from  regions  of  changing  mesh  size.  This  means  that  refining  the 
mesh  in  the  neighborhood  of  points  of  interest  becomes  Impractical  for  high 
frequency  problems  and  that  some  other  method  of  reducing  the  storage  require- 
ments is  needed. 

1.1.2.  PROBLEMS  WITH  NUMERICAL  DAMPING  IN  INTEGRATION  ALGORITHMS 

In  most  structural  dynamics  applications,  the  use  of  Implicit 
unconditionally  stable  algorithms  is  favored.  Such  algorithms  are  the 
Newmark  family  of  operators,  the  Houbolt  method  and  the  Wilson  0 method. 
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These  algorithms  all  provide  numerical  dissipation  to  damp  out  any  previous 
participation  of  higher  mode  response.  The  amount  of  numerical  dissipation 
is  controlled  by  the  time  step.  Reducing  the  time  step  reduces  the  damping 
in  each  mode.  Thus,  to  ensure  accuracy  in  the  higher  modes,  a very  small  time 
step  may  be  needed. 

In  the  Wilson  0 method,  the  only  parameter,  6,  must  be  taken 
greater  than  1.37  to  maintain  unconditional  stability.  It  is  recommended 
(Reference  2)  that  9 be  taken  to  be  1.4  for  which  the  method  is  highly 
dissipative.  To  integrate  a mode  accurately,  it  is  recommended  that  100 
steps  be  used  per  period.  This  would  seem  to  indicate  that  extremely  small 
time  steps  would  be  needed  for  high  modal  accuracy.  Houbolt's  method  is 
generally  considered  to  be  even  more  dissipative  than  the  Wilson  method 
and  the  Newmark  family  is  considered  (Reference  6)  to  be  even  more  dissipa- 
tive than  both  the  Wilson  and  Houbolt  methods.  Thus,  the  Wilson  method  is 
considered  by  many  to  be  the  best  available  unconditionally  stable  algorithm, 
but  is  clearly  unsuited  for  high  frequency  analysis. 

The  small  time  step  and  fine  discretization  needed  for  high  fre- 
quency analysis  by  implicit  analysis  brings  a number  of  attendant  problems. 
These  are  mainly  concerned  with  the  matrix  inversions  associated  with  time 
integration  schemes,  which  will  lead  to  extremely  large  generalized  stiff- 
ness matrices  and  may  make  the  use  of  current  implicit  methods  impractical. 

It  also  suggests  the  possible  use  of  explicit  methods  which  do  not  Involve 
matrix  inversions. 

1.1.3.  PROBLEMS  WITH  MODELING  OF  DAMPING 

Damping  of  structural  vibrations  is  almost  invariably  treated 
according  to  the  viscous  damping  model  in  which  viscous  resisting  forces 
(proportional  to  velocity)  are  assumed  to  act  alongside  the  elastic  restor- 
ing forces  of  the  structure.  One  advantage  of  this  model  is  great  mathe- 
matical simplicity;  the  equations  of  motion  are  very  easy  to  set  up  and, 
under  certain  conditions,  to  solve  provided  that  the  matrix  of  damping  coef- 
ficients is  known.  However,  these  coefficients  are  not  directly  measurable; 
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they  can  only  be  inferred  from  measurements  of  modal  damping  ratios  (since 
it  is  these  ratios  that  deten  ine  the  amplification  factor  at  resonance,  the 
logarithmic  decrement  of  free  vibrations,  etc.).  Experiments  generally  show 
that  the  damping  ratio  in  a homogeneous  structure  varies  little  from  mode  to 
mode.  On  the  other  hand,  the  viscous  damping  model,  with  the  damping  coef- 
ficients regarded  as  material  properties  of  the  st  'uctural  elements  (much  as 
the  elastic  properties),  predicts  that  the  damping  ratio  increases  with 
natural  frequency.  Thus,  the  high  frequency  modes  tend  to  be  damped  out 
(independently  of  the  aforementioned  numerical  damping).  Therefore,  this 
model  is  not  particularly  suited  to  the  study  of  high  frequency  motion. 

A mode-independent  damping  ratio  may  be  assumed  in  the  viscous 
damping  model  only  if  the  mode-superposition  method  is  used  and  this  is 
practicable  only  with  a limited  number  of  degrees  of  freedom.  Other  models 
of  structural  damping  include  the  "hysteretic"  (Reference  7)  and  the 
Rosenblueth-Herrera  (Reference  8)  models,  both  nonlinear,  and  hence  diffi- 
cult to  apply  to  multi  degrce-of-f reedom  systems. 

1.2.  OBJECTIVE 

In  this  section  certain  results  obtained  in  preliminary  investiga- 
tions of  methods  for  the  high  frequency  vibrations  of  structural  systems  are 
described.  To  gain  insight  into  the  problems  which  arise  in  this  area  a series 
of  idealized  models  which  incorporate  the  characteristics  of  an  embedment 
medium,  a structure  and  internal  equipment  was  considered.  The  equipment  has 
natural  frequencies  which  are  high  compared  with  the  fundamental  frequency  of 
the  structure.  The  objective  of  this  analysis  is  to  compare  exact  solutions  to 
the  idealized  model  with  appropriate  solutions  which  may  be  obtained  by 
standard  structural  dynamics  techniques  to  give  an  understanding  of  the  degree 
of  refinement  needed  in  the  modeling  and  the  numerical  technique.  In  the 
course  of  this  study,  certain  results  on  the  vibration  of  equipment  have  been 
obtained  which  appear  to  be  new. 

The  combined  structure-equipment  systems  considered  here  are  such 
that  the  equipment  is  relatively  light  and  has  a natural  frequency  much  above 
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the  fundamental  frequency  of  the  structure.  This  is  In  contrast  to  earlier 
work  of  Geers  and  UeRuntz  (Reference  9),  for  example,  who  consider  struc- 
tural systems,  typical  of  submarines,  where  the  equipment  is  massive,  may  in 
fact  be  heavier  than  the  structure,  and  has  a natural  frequency  lower  than 
the  structural  frequency.  Three  systems  have  been  considered,  two  cf  which 
while  idealized  Incorporate  the  characteristics  of  an  embedment  medium,  a 
structure,  and  attached  equipment.  The  third  model  is  more  representative 
of  silo-like  structures  involving  an  embedded  cylindrical  shell  capped  by  a 
circular  plate  to  which  the  equipment  is  attached. 

Results  presented  here  appear  to  suggest  conclusions  somewhat 
different  from  previous  investigations  such  as,  for  example,  Scavuzzo  (Ref- 
erence 10)  and  Newmark  (Reference  11)  who  have  suggested  that  the  problem  of 
equipment  with  a natural  frequency  close  to  a frequency  of  the  structure  is 
not  a serious  one.  The  present  results  indicate  that  such  a situation,  re- 
ferred to  as  "tuning,"  may  in  fact  be  quite  serious. 

A deficiency  of  the  analysis  for  these  three  models  is  the  neglect 
of  material  damping  (radiation  damping  is  included).  One  objective  is  to 
incorporate  into  the  analysis  damping  models  which  are  realistic  at  high  fre- 
quencies. Many  of  the  difficulties  attendant  on  the  use  of  viscous  damping 
in  structural  analysis  of  Impact  problems  or  high  frequency  analysis  can  be 
alleviated  by  the  use  of  general  viscoelastic  models  to  represent  the  response 
of  the  system.  However,  the  degree  of  complication  which  results  could  be 
considerable  for  a large  system,  and  It  would  be  almost  impractical  to  obtain 
the  necessary  data  to  define  the  response.  What  is  needed  is  some  model 
which  has  the  simplicity  of  the  viscous  model  which  can  be  generally  defined 
by  specifying  a damping  factor  for  each  mode.  In  this  Report  one  example  of 
a model  is  described. 

1.3.  SCOPE 

To  accomplish  these  objectives,  the  Interaction  between  a 
structure-foundation  system  and  an  attached  item  of  equipment  has  been  con- 
sidered through  the  use  of  simple  Idealized  models.  In  the  first  model  a 
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semi-infinite  rod  represents  the  foundation  or  soil,  a finite  rod  represents 
the  structure  and  a spring-mass  system  represents  the  equipment.  The  purpose 
of  the  analysis,  given  in  detail  in  Section  2 is  to  determine  the  acceleration 
amplification  factor  of  the  mass  representing  the  equipment  when  the  system 
is  subject  to  an  enforced  oscillatory  input  of  frequency  oj  at  large  distances 
from  the  foundation-structure  interface.  The  amplification  factor  at  any 
frequency  u)  depends  on  the  following  parameters:  y — the  ratio  of  the  mass  of 
the  equipment  to  that  of  the  structure,  — the  natural  frequency  of  the 
equipment  on  a fixed  support,  0)* — the  fundamental  frequency  of  the  structure 
on  a fixed  foundation  and  a parameter  R which  determines  the  reflection- 
transmission  characteristics  of  the  interface.  A value  of  R less  than  1 
means  that  an  unloading  wave  is  reflected  back  into  the  medium,  greater  than 
1 means  that  a loading  wave  is  reflected. 

To  estimate  its  value  an  unbraced  frame  building  located  on  a 
soil  medium  in  which  a horizontally  polarized  shear  wave  is  propagating 
vertically  has  been  considered.  As  an  example  of  such  a building,  the 
building  discussed  in  Reference  12  on  an  alluvial  soil  was  used  and  a value 
of  R = 0.1  was  obtained.  An  unbraced  frame  is  a relatively  flexible  struc- 
ture, for  a stlffer  building  the  value  of  R would  increase.  When  R = 1,  a 
wave  is  transmitted  without  change  into  the  structure,  the  structure  acting 
in  effect  as  an  extension  of  the  medium. 

The  amplitude  factor  as  a function  of  forcing  frequency  w for  a 
system  with  R = 0.1  was  found  to  have  the  following  characteristics.  For 
small  values  of  the  mass  ratio,  y,  the  resonances  of  the  system  below  the 
natural  frequency  of  the  equipment  are  shifted  below  the  natural  frequencies 
of  the  structure  alone  and  those  above  are  shifted  above  the  structure  fre- 
quencies. If  the  equipment  frequency  is  not  close  to  a structure  frequency, 
an  additional  resonance  is  interpolated  between  those  structural  frequencies 
immediately  above  and  below  the  equipment  frequency.  The  results  for  an 
untuned  system  are  shown  in  Figures  5a  and  5b.  In  these  figures  the  amplifica- 
tion factor  is  shown  for  a case  where  R = 0.1  and  the  natural  frequency  of 
the  equipment  lies  between  4 and  5 times  the  fundamental  frequency  of  the 
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When  the  equipment  resonance  falls  on  a structure  resonance  (i.e., 
the  system  Is  tuned),  a double  peak  occurs  with  the  two  maxima  close  together. 

It  is  easy  to  show  that  the  distance  between  the  two  peaks  is  proportional  to 
/y.  An  interesting  point  is  that  if  the  equipment  is  not  tuned  to  the  struc- 
ture the  dominant  peak  is  virtually  at  the  equipment  resonance;  however,  if 
it  is  tuned,  the  value  of  the  amplification  factor  at  the  equipment  frequency 
is  reduced  and  is  proportional  to  y.  The  peaks  on  either  side,  however,  are 
now  important.  These  results  are  shown  in  Figure  6a  and  also  Figure  6b  as  an 

expanded  curve  near  the  equipment  resonance  when  the  equipment  mass  ratio  is  i 

very  small.  | j 

i 

This  double  peak  corresponds  to  the  occurrence  of  two  simple  poles  i 

I 

of  the  system  transfer  function  in  the  transform  domain  being  very  close  to-  j 

gether  and,  in  the  case  of  very  light  mass,  to  essentially  a double  pole.  I 

The  implications  of  this  double  peak  on  the  transient  response  of  the  system  '\ 

could  be  considerable.  It  is  easy  to  show  that  the  presence  of  a double  pole  i 

leads  to  transient  build-up  which  is  proportional  to  t in  a condition  of  i 

1 

resonance.  The  tuned  situation  thus  appears  to  be  more  critical  even  than  J 

normal  resonance,  in  contrast  to  earlier  results.  For  example,  Scavuzzo  ] 

(Reference  10)  using  somewhat  similar  modeling  involving  an  infinite  bar  to  ‘ 

1 

which  a cantilevered  mass  was  attached,  found  by  numerical  techniques  that 
the  response  at  the  frequency  of  the  vibrating  mass  is  less  when  the  mass  is 

tuned  than  when  it  is  not  and  concluded  that  tuning  is,  therefore,  not  a j 

critical  situation.  Herein  the  identical  result  is  obtained  but  the  inter-  I 

pretation  is  quite  different.  The  response  at  the  frequency  of  the  equipment  j 

when  the  system  is  not  tuned  is  very  large  and  when  tuned,  is  much  less.  I 

This  is  clearly  indicated  in  Figures  5 and  6.  However,  the  fact  is  that  a . 

numerical  method  might  well  overlook  the  very  large  amplification  factors 
on  each  side  of  the  tuned  equipment  frequency  when  the  system  is  tuned. 

These  twin  peaks  may  lead  to  amplification  over  a range  of  frequencies  on 
either  side  of  the  equipment  frequency,  thus  pre:  anting  a large  target  to 
resonant  conditions. 


The  results  of  this  analysis  on  the  elementary  models  may  be 
summarized  as  follows: 


a.  In  the  structure-equipment  case,  the  frequencies  of  the  system 
are  shifted  such  that  those  below  the  equipment  frequency  are  lower,  those 
above  higher. 

b.  As  the  mass  ratio  becomes  small,  the  frequencies  tend  to 
the  structural  frequencies  and  the  equipment  frequency. 

c.  If  the  equipment  frequency  is  not  tuned  to  the  structural 
frequencies,  the  response  curve  for  light  mass  is  roughly  the  superposition 
of  the  structural  response  and  the  equipment  response  curves  with  little 
interaction. 

d.  If  the  equipment  frequency  is  tuned  to  a structural  frequency, 
there  are  two  frequencies  close  together  on  either  side  of  the  tuning  fre- 
quency. The  response  magnitude  Intermediate  between  these  is  proportional  to 
1/Y  while  the  spread  between  the  frequencies  is  proportional  to  /y.  There 

is  thus  a spread  of  high  amplification  over  that  band. 

e.  In  the  transient  problem  the  double  peak  resonance  is  equiva- 
lent to  a multiple  pole  in  the  limit  as  Y "*■  0.  The  resonant  response  around 
the  tuning  frequency  will  grow  like  t^  rather  than  t.  Thus,  in  the  undamped 
case  this  will  dominate. 

f.  Since  the  spread  between  the  peak  frequencies  in  the  tuned 
case  is  proportional  to  the  square  root  of  the  mass  ratio,  the  effect  is  a 
band  of  resonant  frequencies  which  vanishes  slowly  with  decreasing  mass 
ratio  and,  thus,  for  light  masses  offers  a substantial  target  for  sympathetic 
oscillation. 

A third  model  more  representative  of  silo-like  structures  was 
analyzed  and  it  was  found  that  essentially  the  same  type  of  response  is 
developed  in  more  realistic  models  of  equipment-structure  systems  and  that 
the  phenomenon  of  the  double  peak  in  a tuned  system  is  not  confined  to  the 
elementary  model.  The  structure  considered  here  is  illustrated  in  Figure  7. 

It  comprises  a long  cylindrical  shell  topped  by  a flat  plate.  To  maintain 
the  analytical  simplicity  afforded  by  axial  symmetry  in  the  response  of  the 
plate  and  cylinder  the  equipment  is  idealized  as  a spring-supported  mass 
connected  to  the  structure  at  the  center  of  the  plate. 
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The  system  Is  assumed  to  be  excited  by:  (a)  sinusoidal  displace- 
ment input  at  the  base  of  the  shell  which  is  far  removed  from  the  plate;  or 
(b)  sinusoidally  varying  pressure  loading  on  the  plate;  or  (c)  an  oscillatory 
point  load  applied  to  the  center  of  the  plate.  These  excitations  maintain 
axial  symmetry  and  simulate  ground  shock,  air  shock  and  shaker  excitations, 
respectively.  Other  cases  which  do  not  possess  axial  symmetry  could  be  treated 
but  with  a substantial  increase  in  analytical  complexity. 

The  algebraic  complexity  of  the  analysis  is  considerably  greater 
than  for  the  elementary  models  studied  in  the  earlier  report  and  the  fact 
that  the  plate  is  governed  by  a fourth  order  equation  while  the  bar  equation 
is  second  order  gives  rise  to  a number  of  differences  in  the  results.  How- 
ever, the  essential  features  are  the  same.  The  double  resonance  occurs  as 
before  when  the  system  is  tuned.  The  double  peak  resonance  is  as  before 
dominant  in  the  amplification  diagram  and  it  is  found  that  as  the  equipment 
mass  becomes  small  the  spread  between  the  peaks  decreases  at  a rate  propor- 
tional to  the  square  root  of  the  mass  ratio.  Results  derived  for  the  spread 
between  the  peaks  as  a function  of  the  structure  frequency  to  which  the  equip- 
ment is  tuned  indicate  that  the  phenomenon  is  more  important  at  frequencies 
above  the  fundamental  and  the  effect  appears  to  decrease  at  very  high  fre- 
quencies of  tuning.  The  phenomenon  is  thus  an  intermediate  frequency  one  in 
this  case,  but  it  is  clear  that  it  may  arise  in  a variety  of  structure- 
equipment  systems. 

To  accomplish  the  objective  of  incorporating  simple,  but  realistic 
damping,  a model  based  on  the  simple  linear  viscoelastic  material  known  as 
the  Maxwell  material  has  been  developed.  The  conventional  viscous  model  of 
structural  damping  is  analogous  to  the  Kelvin  model  of  viscoelasticity  and 
the  other  simple  model  is  that  denoted  as  the  Maxwell  model.  It  is  well  known 
that  the  Kelvin  model  is  inadequate  for  wave  propagation  problems  since  the 
model  cannot  transmit  waves;  it  implies  an  Infinite  wave  speed  and  leads  to 
parabolic  equations  of  motion.  The  Maxwell  model  on  the  other  hand  leads  to 
hyperbolic  equations  of  motion  and  transmits  waves  with  the  elastic  wave 
speed  and  some  attenuation.  Thus,  it  seems  likely  that  the  Maxwell  model 


15 


when  converted  to  the  structural  problem  might  be  useful  for  impact  problems 
and  high  frequency  analysis.  It  should  be  noted  that  the  same  number  of  para- 
meters are  needed  to  define  the  Maxwell  model  as  are  needed  for  the  Kelvin 
model  and  that  the  use  of  modal  damping  factors  would  also  suffice. 

The  most  effective  way  to  illustrate  the  behavior  of  the  damping 
model  suggested  here  is  to  study  the  response  of  a single  degree-of-f reedom 
system.  This  is  done  in  detail  in  Section  4 where  the  steady  state  response 
of  a one  degree-of-f reedom  system  under  sinusoidal  loading  is  determined. 

The  steady  state  dynamic  load  factor,  that  is  displacement  per  unit  load,  is 
shown  in  Figure  14.  In  the  region  of  resonance  the  response  is  similar  to 
that  for  Kelvin  damping  but  near  zero  frequency  the  result  grows  without 
bound  as  the  frequency  tends  to  zero.  This  is  to  be  expected  since  the 
Maxwell  model  behaves  as  a fluid  unaer  static  loading. 

For  a multi  degree-of-f reedom  system,  the  equations  of  motion  for 
the  Maxwell  model  can  be  put  into  the  same  form  as  those  for  conventional 
viscous  damping  with  a modified  forcing  function  vector.  This  Implies  that 
a computer  program  designed  for  the  conventional  dynamic  analysis  of  struc- 
tural systems  can  be  employed  with  slight  modification  for  the  same  type  of 
problem  but  with  the  damping  represented  by  the  Maxwell  viscoelastic  model. 

A significant  difference  that  exists  between  the  two  models  is  in 
the  initial  conditions  which  need  to  be  specified.  The  conventional  viscous 
model  requires  that  the  displacement  and  velocity  at  the  initial  time  be 
known  and  these  quantities  completely  determine  the  initial  internal  forces. 
With  a Maxwell  model,  however,  these  quantities  do  not  determine  the  initial 
value  of  the  Internal  force.  In  fact,  this  must  be  specified  in  addition  to 
the  displacement  and  the  velocity.  This  aspect  of  the  Maxwell  model  suggests 
that  it  might  indeed  be  more  suitable  than  the  Kelvin  model  for  impact  prob- 
lems and  similarly  for  high  frequency  analysis.  While  a structural  system  at 
rest  acquires  an  Initial  velocity  following  an  Impact,  it  is  not  reasonable 
to  assume  that  this  also  results  in  non-zero  values  for  resisting  forces. 

The  use  of  a Kelvin  model  does  Incorporate  such  an  assumption  and  leads  to 
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high  accelerations  and  low  displacements  In  the  higher  modes.  In  contrast, 
the  Maxwell  model  results  in  more  flexible  response  and  gives  lower  accelera- 
tions and  higher  displacements.  With  mass  proportional  damping,  modal  damp- 
ing is  inversely  proportional  to  modal  frequency  and  both  models  give  essen- 


l tially  the  same  answers  for  the  low  damping  ratios  resulting  from  sucli  an 

i assumption. 

I 


SECTION  2 


ANALYSIS  OF  EQUIPMENT-STRUCTURE  INTERACTION  IN  IDEALIZED  SYSTEMS 
2.1.  FIRST  ANALYTICAL  MODEL 

As  an  approach  to  a qualitative  understanding  of  the  problem  of 
equipment  response  to  ground  shock,  the  following  elementary  model  is  con- 
sidered. The  structure  is  represented  by  a uniform  rod  of  area  a,  Young's  I 

modulus  E,  density  p of  semi-infinite  extent  -“  < x < 0 to  which  a mass  m is  i 

attached  by  a spring  k which  represents  the  equipment.  Displacement  of  the 

rod  is  given  by  u (x,  t),  -«>  < x < 0,  and  of  the  equipment  by  w (t).  The  | ’ 

rod  is  subject  to  a forced  motion  u at  distances  remote  from  x = 0.  ■ 

° i 

Figure  1 Illustrates  the  system.  ^ 

I 

] 

The  equations  of  motion  are  j 

m w = - F where  F = k(w-u(o,t)),  '''■ 

1 

I 

so  that  m w + k w = k u(o,t),  (2-1)  > 

and  p u,  ^ = E u,  with  F = E a u,  (o,t).  (2-2)  | 

tt  XX  X ] 

i 

The  solution  in  the  rod  is  given  by 

u(x,t)  = c^  = E/p  (2-3)  i 

i 

the  first  term  representing  the  ground  shock  input  and  the  second  the  back  j 

] 

reflected  wave  from  the  equipment  structure  Interface.  It  follows  that 

F = E a la)/c  (A-B)e^“‘^ 
and  taking  w = results  in 

(-iwi)  + k)We  = k (A+B)e 

and  E a iw/c  (A-B)e^“‘^  = k {A+B-W}  e^^*^  I 

J 

I 

Thus,  there  are  two  equations 

k j 

w = (A+B)  (2-4) 

K— mu)  j 


$ 
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and 


W = A+B  - (A-B) 


(2-5) 


relating  W,  A,  B.  A is  assumed  known  and  W to  be  determined. 
Defining 


D = l/(l-a)^/a)  ) with  = k/m 
o o 


it  follows  that 


W = - 2i 


Eaoj/kc 


(1-1/D)  - i Eato/kcD 


(2-6) 


The  comparison  of  the  maximum  acceleration  of  the  equipment,  i.e,  u)^  W to  the 
maximum  acceleration  induced  by  the  ground  shock  A is  given  by 

2a(i)/a) 


{ (1-1/D) ^ 

o 2 


= 2 a 


(2-7) 


where  (2  = oi/u)  and  a = pac/mw  . 

o o 

(a  is  a sort  of  impedance  ratio,  or  could  be  Interpreted  as  a mass  ratio  if 

c/w  is  taken  as  a length, 
o 

The  function  (2^  + a^(l-(2^)^  has  a minimum  where 
= (2a^-l)/2a^  if  > 1/2 
and  the  value  of  this  minimum  is 

(4a^-l)/4a^ 

Thus, 


= 2,  < 1/2 


max 


4aV(4a^-l)^'^^  if  > 1/2 


(2-8) 

This  result  is  shown  in  Figure  2 in  which  |w/A|  is  plotted  as  a function  of 
J2  for  various  values  of  a. 
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The  conclusion  is  that  the  amplification  factor  has  its  maximum 
value  at  a frequency  shifted  below  the  equipment  frequency  and  that  amplifica- 
tion tends  to  a maximum  value  of  2a  for  large  a.  Large  a corresponds  to 
light  equipment;  thus,  amplification  can  be  very  large  for  light  equipment. 

The  disadvantages  of  this  model  are  that  the  structure  is  not  clearly  defined. 

It  is  interesting  to  note  that  this  result  has  an  analogy  with  the 
single  degree-of-f reedom  damped  oscillator.  The  mass  ratio  a in  this  analysis 
corresponds  to  1/(2?)  where  f is  the  fraction  of  critical  damping  in  the  damped 
oscillator.  The  amplification  factors  in  the  two  cases  differ  by  a factor  of 
2,  which  arises  in  the  present  case  due  to  the  doubling  (caused  by  wave  reflec- 
tions) at  the  right  end  of  the  bar  of  the  motion  imposed  at  the  left  end  of 
the  bar.  Critical  damping,  \ = 1,  corresponds  here  to  a mass  ratio  a = 1/2, 
at  which  value  the  free  vibration  response  becomes  non-osclllatory . The 
amplification  factor  at  any  frequency  ratio  12  is  less  tlian  that  at  zero  fre- 
quency for  values  u < 1//2  even  though  the  free  vibration  response  is  oscil- 
latory wtien  a > 1/2.  It  is  also  useful  to  note  that  the  amplification  factor 
at  12  = 1,  i.e.,  at  resonance  of  the  undamped  system,  is  less  than  the  zero 
frequency  value  when  a < 1,  which  corresponds  in  the  analogy  to  a damping 
ratio  of  5 = 1/2.  Thus,  in  resonance  testing  of  buried  structures  where 
Interest  is  concentrated  on  the  frequency  ratio  near  the  resonance  of  the  un- 
buried structure,  the  resonance  of  the  buried  structure  would  appear  to 
vanish  for  an  effective  damping  ratio  of  around  1/2,  equivalent  to  a mass 
ratio  of  around  1.0,  even  though  this  buried  system  has  an  oscillatory 
response  in  free  vibration.  This  may  be  an  explanation  for  the  unanticipated 
disappearance  of  resonance  peaks  in  resonance  testing  of  burled  structures. 

2.2.  ADVANCED  ANALYTICAL  MODEL 

A disadvantage  of  the  previous  model  is  that  the  mass  of  the  struc- 
ture is  not  easily  specified.  In  order  to  obtain  a definite  structure  mass, 
a model  in  which  a semi-infinite  rod < x < 0,  represents  the  embedment 
medium,  a finite  rod  0 < x < £ represents  the  structure  and  a spring  mass 
system  represents  the  equipment  as  shown  In  Figure  3 is  considered.  The  dis- 
placement of  the  semi-infinite  rod  will  be  denoted  by  u^(x,t),  -<»  < x < 0,  that 
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of  the  finite  rod  by  u,(x,t),  0 < x < ?..  As  before  the  displacement  of  the 
mass  will  be  denoted  by  w(x,t),  and  it  is  assumed  that  the  system  is  subject 
to  an  enforced  input  at  large  distances  from  the  interface.  Thus,  the  solu- 
tion involves  three  parts: 

u^  = ^ ^ia)(t+x/c^)  ^2-9) 


U2  = (A2  sin  tjx/c2  + 82  cos  Uix/c^)  e 


(2-10) 


m w + k w = k U2(£,t),  where  c^  = = E2/P2  (2-11) 


with  the  conditions 


(2-12) 


(2-13) 


Thus , 


Uj^(0,t)  = U2(0,t); 


k (w-u2(i.,t))  = £232^2 


Ai  + Bi  = B2 


E2aj^(-ia)/Cj^Aj^  + iw/Cj^B^)  = 


kW-k  (A^  sin  u^Z/c^  + cos  taZ/c^)  = E2a2{w/c2A2  cos  iiiZ/c^-ui/ c^B^  sin  u^Z/c^} 

As  before,  A^  is  assumed  specified  and  W is  solved  for  by  elimination  of  B^^, 
A2,  and  B2.  After  much  manipulation,  it  is  found  that: 


where 


sin  tiiZ/cj  + r cos  oiJc/c- 
= 2d A^ 


(D-l)  sin  UiZ/c^  - E2a2w/kc2  cos 
^ (1-D)  cos  iiiZ/c^  - E2a2to/kc2  sin  iiiZ/c^ 


R = E2a2C^/Ej^aj^C2  • 


(2-14) 


(2-15) 


(2-16) 


The  amplification  factor,  namely,  the  absolute  magnitude  of  W/A^,  is  given  by 

u ^0^9*^  E a 0) 

— = 2 — {((D-1)  sin  wJl/c^ ■ cos 


+ ((1-D)  cos  w£/c„ sin  tu^/c.)^} 


(2-17) 


Noting  that  the  total  mass  of  the  structure  M is  given  by  P2^2^  that  the 
fundamental  frequency  of  the  structure  w*  is  '^c^/1%  a mass  ratio  y = '"/M  is 
introduced  finally  yielding  the  expression 


where 


Y = 2 {P^  + 


T,  /i2/2x  TTW  TTCO  .TTW 

P = (1-0)  /(jo  ) cos  ^ - Y ^ sin  ^ —r 

o 2 0)*  2 0)*  2 a)* 


_ 2/2\  . TTtO.^Tia)  DO) 

Q = (1-0)  /o)  ) sin  -Jr  — X + Y TT  ^ cos  — — 
o 2 0)*  2 0)*  2 0)* 


(2-18) 


(2-19) 


(2-20) 


In  this  equation  the  parameters  Y.  w*  can  be  easily  given 
numerical  values,  but  the  parameter  R is  more  difficult  to  determine.  To 
estimate  its  value  the  example  is  considered  of  a shear  building  with  column 
stiffness  k,  story  masses  m,  story  height  h and  story  module  b x d,  located 
on  a soil  medium  with  density  and  shear  modulus  in  which  a horizontally 
polarized  shear  wave  is  propagating  vertically. 


The  equation  of  motion  for  the  i mass  of  a module  of  the  shear 


building  is 


m = - 2k  u^  + k + k 


which  leads  to  a continuum  equation  of  the  form 

m u, = k h^  u, 

tt  XX 

Dividing  both  sides  by  b d h gives  an  equivalent  density  p^  = 
m/bdh  and  equivalent  modulus  = 2kh/bd  and  thus  in  the  parameter  R = 

E a.,c,/E  a,c_,  E, , c,  is  replaced  by  the  shear  modulus  G and  shear  speed 


I 
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1/2  1/2 
c = (G  /p  ) of  the  soil  and  E.,c,  by  E,  and  c,  = (E, /p,  ) and  take 
sss  ^2  b bbb 

a,  = a^.  Thus,  R = kh/bd  (m/kh^) '^^/G  c . An  alternative  form  is  to  write 
12  s s 

R in  the  form  p,  c. /p  c , and  utilize  information  on  the  fundamental  period 
D b s s 

of  the  structure.  For  example,  a shear  building  (Reference  12)  of  total 
height  230  feet  has  a period  of  2.2  seconds.  The  wave  speed  is  thus  AOO  feet/ 
second.  The  equivalent  density  of  this  building  is  estimated  to  be  50  pounds/ 
feet^  based  on  a story  weight  of  250,000  pounds,  a story  height  of  12  feet  and 
a 20  feet  by  20  feet  module.  Considering  a soil  density  of  100  pounds/feet^ 
and  a soil  shear  wave  speed  of  2,000  feet/second,  R = 0.1. 

The  number  R determines  the  reflection  transmission  characteristics 
of  the  interface.  A value  of  R less  than  1 means  that  an  unloading  wave  is 
reflected  back  into  the  medium,  greater  than  1 means  a loading  wave  is  reflected. 
A shear  or  unbraced  frame  building  is  a flexible  structure;  for  a stiffer  build- 
ing the  value  of  R would  increase.  When  R = 1,  a wave  is  transmitted  without 
change  into  the  structure,  the  structure  acting  in  effect  as  an  extension  of 
the  medium.  Thus,  in  the  case  of  R = 1,  the  result  should  reduce  to  that  of 
the  earlier  model,  and  it  is  easy  to  see  that  it,  in  fact,  does.  It  is,  there- 
fore, useful  to  consider  the  case  of  small  R since  for  R near  1 the  previous 
results  will  hold. 

For  small  values  of  R the  peaks  of  the  amplification  function  |w/A^| 
will  be  close  to  the  zeroes  of  the  term  P.  These  are  given  by 

2 2 ^ 

When  m = 0,  the  zeroes  are  those  of  cos  ^ namely  (2n-l)/2  it,  which  cor- 

l (jj* 

respond  to  a bar  fixed  at  one  end  and  free  at  the  other.  The  zeroes  for  finite 
values  of  m are  shifted  lower  for  those  below  cj  = (the  natural  frequency  of 
the  equipment)  and  shifted  higher  for  those  above. 

When  X tan  x and  l/y  (l-3^x^)  where  3 = 2u)*/7ra)^  are  plotted  it 
is  seen  that  the  zeroes  of  P are  the  intersections  of  these  curves.  For  small 
values  of  m,  it  is  clear  from  Figure  A that  for  the  untuned  case  of  the  three 
zeroes  nearest  the  equipment  resonance  one  is  slightly  below  the  structure 
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resonance  immediately  below  the  equipment  resonance,  one  is  close  to  the 
equipment  resonance  and  one  is  slightly  above  the  structure  resonance  im- 
mediately above  the  equipment  resonance.  When  the  equipment  resonance  falls 
on  a structure  resonance  (i.e.,  the  system  is  tuned),  a double  root  occurs 
with  two  intersections  very  close  together.  It  is  easy  to  show  that  the 
distance  between  the  two  roots  is  proportional  to  /y.  An  interesting  point 
is  that  if  the  equipment  is  not  tuned  to  the  structure  the  dominant  peak  is 
virtually  at  the  equipment  resonance;  however,  if  it  is  tuned,  the  value  of 
the  amplification  factor  at  o)  = is  reduced  and  is  proportional  to  l/y* 

The  peaks  on  either  side  are  now  important.  These  results  are  shown  in 
Figures  5,  6a  and  6b.  In  Figure  5 the  result  is  shown  when  the  system  is 
not  tuned  and  in  Figure  6a  the  result  when  it  is.  Figure  6b  shows  an  expanded 
curve  near  the  equipment  resonance  when  the  equipment  mass  ratio  is  very 
small. 

The  implications  of  this  double  root  on  the  transient  response 
of  the  system  are  considerable.  It  is  easy  to  show  that  the  presence  of  a 
double  root  leads  to  transient  buildup  which  is  proportional  to  t^  in  a condi- 
tion of  resonance.  The  tuned  situation  thus  appears  to  be  more  critical 
even  than  normal  resonance  in  contrast  to  earlier  results.  For  example, 
Scavuzzo  (Reference  10),  using  a somewhat  similar  modeling  involving  an 
infinite  bar  to  which  a cantilevered  mass  is  attached,  found  by  a numerical 
technique  that  the  response  at  the  frequency  of  the  vibrating  mass  is  less 
when  the  mass  is  tuned  than  when  it  is  not  and  concluded  that  tuning  is, 
therefore,  not  a critical  situation.  The  result  obtained  here  is  the  same; 
the  response  at  the  frequency  of  the  equipment  when  the  system  is  not  tuned 
is  very  large  and  when  tuned  is  much  less.  This  is  clearly  indicated  in 
Figures  5 and  6.  However,  the  fact  is  that  a numerical  method  might  well 
overlook  the  very  large  amplification  factors  on  each  side  of  the  tuned 
equipment  frequency  when  the  system  is  tuned.  These  twin  peaks  may  lead  to 
amplification  over  a range  of  frequencies  on  either  side  of  the  equipment 
frequency,  thus  presenting  a larger  target  to  resonant  conditions. 


SECTION  3 


L 


ANALYSIS  OF  EQUIPMENT-STRUCTURE  INTERACTION  IN  SILO-LIKE  STRUCTURES 

This  section  describes  the  application  of  the  techniques  developed 
in  the  previous  section  on  high  frequency  structure-equipment  interaction  to 
the  vibration  of  an  equipment-structure  system  which  is  more  representative 
of  silo-like  protective  structures.  The  purpose  of  this  is  to  show  that 
essentially  the  same  type  of  response  is  developed  in  more  realistic  models 
of  equipment-structure  systems  and  that  the  phenomenon  of  the  double  peak  in 
a tuned  system  is  not  confined  to  the  elementary  model.  The  structure  con- 
sidered here  is  Illustrated  in  Figure  7.  It  comprises  a long  cylindrical 
shell  topped  by  a flat  plate.  To  maintain  the  analytical  simplicity  afforded 
by  axial  symmetry  in  the  response  of  the  plate  and  cylinder  the  equipment  is 
taken  to  be  a spring  supported  mass  connected  to  the  structure  at  the  center 
of  the  plate. 

The  system  is  taken  to  be  excited  by:  (a)  sinusoidal  displace- 
ment input  at  the  base  of  the  shell  which  is  far  removed  from  the  plate;  or 
(b)  sinusoidally  varying  pressure  loading  on  the  plate;  or  (c)  an  oscillatory 
point  load  applied  to  the  center  of  the  plate.  These  excitations  maintain 
axial  symmetry  and  simulate  ground  shock,  air  shock  and  shaker  excitations, 
respectively.  Other  cases  which  do  not  possess  axial  symmetry  could  be 
treated,  but  with  a substantial  Increase  in  analytical  complexity. 

The  algebraic  complexity  of  the  analysis  is  considerably  greater 
than  for  the  elementary  models  studied  in  the  earlier  report  and  the  fact 
that  the  plate  is  governed  by  a fourth  order  equation  while  the  bar  equation 
is  second  order  gives  rise  to  a number  of  differences  in  the  results.  How- 
ever, the  essential  features  are  the  same.  The  double  resonance  occurs  as 
before  when  the  system  is  tuned.  The  double  peak  resonance  is  as  before 
dominant  in  the  amplification  diagram  and  it  is  found  that  as  the  equipment 
mass  becomes  small  the  spread  between  the  peaks  decreases  at  a rate  propor- 
tional to  the  square  root  of  the  mass  ratio.  Results  have  also  been  derived 
between  the  peaks  as  a function  of  the  structure  frequency  to  which  tlie 
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equipment  is  tuned.  These  results  indicate  that  the  phenomenon  is  more 
important  at  frequencies  above  the  fundamental  and  the  effect  appears  to 
decrease  at  very  high  frequencies  of  tuning.  The  phenomenon  is  thus  an 
intermediate  frequency  one  in  this  case,  but  it  is  clear  that  it  may  arise 
in  a variety  of  structure-equipment  systems. 

3.1.  GOVERNING  EQUATIONS  OF  EQUIPMENT-STRUCTURE  SYSTEM 

The  system  to  be  considered  consists  of  three  elements,  the 
semi-infinite  cylindrical  shell,  the  circular  plate  and  the  spring  supported 
mass  which  simulates  the  equipment. 

The  plate  equations  used  here  are  those  of  elementary  Kirchhoff 
plate  theory;  all  quantities  are  assumed  to  be  functions  of  the  radius  r 
and  time  t.  The  surface  tractions  on  the  plate  are  positive  in  the  direction 
of  positive  transverse  displacements  w of  the  middle  surface  of  the  plate. 

The  transverse  shear  force  is  denoted  by  Q^  and  is  given  in  terms  of  the 
displacements  by 

Q = - D (V^w),  (3-1) 

r p r 

where  D is  the  plate  stiffness. 

P 

The  equation  of  motion  for  the  plate  is 

D V**  w (r,t)  + y w (r,t)  •=  p (r,t),  r > 0 (3-2) 

P P 

where  p (r,t)  is  the  externally  applied  distributed  load  on  the  plate  which 
is  intended  to  approximate  alrblast  effects  on  the  surface  of  the  plate. 

The  shell  equations  to  be  used  are  based  on  axially  symmetric 
membrane  behavior.  In  addition,  the  influence  of  lateral  inertia  is  neglected. 
The  shell  is  also  considered  to  be  very  rigid  in  bending  so  that  only  its 
longitudinal  behavior  is  of  consequence.  This  also  implies  that  the  plate 
edge  is  constrained  against  rotation,  i.e.,  one  appropriate  plate  boundary 
condition  is 

w,^  (a,t)  = 0.  (3-3) 


Oil 


In  order  to  simulate  the  effect  of  ground  shock  the  longitudinal 
displacement  of  the  shell  Is  given  by 

u = A + B (3-4) 

s s s 

where  c = /e  /p  . The  first  term  represents  the  prescribed  ground  shock 
s s s 

input  at  the  shell  base  and  the  second  the  back  reflected  disturbance  at  the 
plate-shell  interface.  When  the  problem  relates  to  air  shock  or  shaker 
excitation  on  the  plate  the  first  term  is  omitted. 

The  equipment  is  taken  to  be  a single  degree-of-f reedom  undamped 
oscillator  with  natural  frequency  = A/m.  The  force  between  the  equipment 
and  the  plate  is  denoted  by  f(t). 


The  interaction  equations  between  the  various  elements  are  as 


follows : 


(a)  The  shell-plate  interaction  is  governed  by  continuity  of 
displacement  and  a balance  between  the  plate  edge  shear  and  the  longitudinal 
force  in  the  shell,  leading  to 


w (a,t)  = u^  (0,t) 


Qr  = ^s  “s  ‘'s 


(b)  The  equipment-plate  interaction  is  governed  by 
f (t)  = k (u-w(0,t)) 


mil  = - f (t)  . (3-8) 

The  appropriate  regularity  and  force  balance  conditions  at  the 
plate  center  are 

w (0,t)  < “ (3-9) 


f 


and 


lim  271  rQ  = - (f+P(t)) 
r ->  0 


(3-10) 


where  P(t)  is  the  applied  shaker  force  at  the  center  of  the  plate. 

Thus,  the  three  possible  excitations  of  interest  are  represented  by 

iwt 


p (r,t)  = p e 

iwt 


(3-11) 

P (t)  = P e"'"'-  (3-12) 

and  exp  [i  (t+x/c^)]]  and  in  the  subsequent  analysis  the  quantities  P,  P, 
and  u)  will  be  considered  as  specified.  The  acceleration  of  the  equipment 
mass  will  be  the  result  of  interest.  Since  the  inputs  are  oscillatory  with 
specified  frequency  o)  it  is  convenient  to  introduce  the  variables  W,  U,  F, 
where 


"!  cot 

w (r,t)  = We  , W=W  (r) 

TT 

u = U e 

iwt 


(3-13) 

(3-14) 


f = F e"”"  (3-15) 

in  terms  of  which  the  governing  equations  of  the  plate  and  the  interactions  are 

D V^W-u  w^W=p,  r>0  (3-16) 


W,^  (a)  = 0 


lim  W < °° 
r -*■  0 


lim  27tr  D (V^W)  , = F + P 

r - 0 P 


W (a)  = A + B 
s s 

D (V^W(a)),  = — (A  -B  ) E t 

p r c s s s s 


(3-17a) 

(3-17b) 

(3-17c) 

(3-17d) 

(3-17e) 


with 


F = mw^  U = k (U-W(O)) 


(3-18) 
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3.2.  A.NALYSIS 


F 


The  general  solution  of  Equation  (3-16)  is 

W (r)  = A J (otr)  + B Y (ar)  + C I (ar) 
o o o 


+ D K (ar)  - 7-^-2 


(3-19) 


where  for  simplicity  we  have  taken  p(r)  = const.  = p and  where 


o'*  = p (U^/D  . 
P P 


(3-20) 


In  this  equation  A,  B,  C and  D are  constants  to  be  determined 

from  the  boundary  conditions  and  J , Y are  Bessel's  functions  of  zero  order 

o o 

of  the  first  and  second  kind,  respectively,  and  I , K are  modified  Bessel 

o o 

functions  of  zero  order.  The  functions  Y (x)  and  K (x)  exhibit  singular 

00 

behavior  at  x = 0 having  representations  near  x = 0 of  the  forms 


Y (x)  = — In  X 
O IT 


(x)  = - In  X,  X = 0. 


Thus,  the  regularity  condition  (17c)  requires 


B - = D 

TT 


leading  to 


W (r)  = A J (ar)  + C I (ar)  + B (Y^(ar)  + - (ar)) *-^2 

^P 

Using  Equation  (3-1)  leads  to 

11m  2TTrQ  = - 11m  2it  D a^r  (A  J,(ar)) 
r ->  0 r - 0 P ^ 


+ C (ar)  + B (Y^(ar)  - - K^(ar)) 


(3-21) 


The  functions  and  exhibit  singular  behavior  at  r = 0 of  the  form 


Y,  (oir)  = - 
1 TOr 


Ki  <ar)  ■ Jf.  r =.  0 


whereas 


(0)  = 0 


and 


(0)  = 0. 


This  Equation  (3-21)  leads  to 


„ F + P 
B = - 2 


8D 

P 


and 


W (r)  = A (ar)  + C (ar)  - ^^2  (Y^  (ar)  + f (ar)) 

P 


p 0) 
P 


When  applied  to  the  above  equation  boundary  condition  (3-17a)  leads  to 


- A J,  + C I,  + ^^-^2  (Y,  + I KJ  = 0 


11  8D  '1  TT  1' 
P 


(3-22) 


where,  above,  and  in  what  follows,  the  arguments  of  the  Bessel  functions  are 
aa. 


The  interaction  Equation  (3-17d)  gives 


A J + C I v,,~2  (Y  + - K ) ^2  = A + B 

o o 8D  a o ir  o p u)  s s 

P P 


(3-2  3) 


and  Equation  (3-17e) 
D 

P 


-D  Ca  Jj  t C - Ijtl,  (Yj  - ^ K^)]  . f (3-24) 

P s 


1(1) 
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Equation  (3-18)  provides  finally 


(-mo)^  + k)U  = k(A  + C-  -^2)  (3-25) 

p to 
P 

Substitution  of  mto^U  for  F in  the  above  system  of  Equations  (3-22  through  3-25) 
leads  to  four  equations  relating  the  four  unknowns  A,  C,  U and  in  terms  the 
input  quantities  A^,  p,  P and  the  frequency  to.  The  important  quantity  is  the 
absolute  acceleration  of  the  mass  |to^u|  and  is  not  necessary  to  obtain  A,  C 
and  By  manipulation  of  Equation  (3-23),  Equation  (3-24)  B^  can  be 

eliminated  obtaining 

A (J  + i6  J,)  + C (I  + ig  IJ  - 
o 1 o 1 


F + P 
8D 


P 


[Y  + ig  Y-  + |-  (K  - ig  K )]  - -^2  = 2 A 

“-O  iTTo  I-*  UtO  S 

P 


(3-26) 


where  g = D ot^  c /toE  t . 

p s s s 

Equations  (3-22)  and  (3-25)  can  now  be  used  to  determine  expres- 
sions for  A and  C in  terms  of  U,  p and  P which  are  substituted  into  Equation 
(3-26)  above  to  determine  a single  expression  for  U in  terms  of  p,  P and  A^. 
The  result  may  be  expressed  in  the  form 


{N  + iR  (-  Q}  U = 2 A + -2-2  {S  + iR  (-  )^^^  T} 

Wj.  s p 0)^  0)^ 

p 


where 


{[j  I,  + I Jj  [1  - - 2] 

•-ol  Cl)-* 


* 1'  S.  -"1  * f ''l>  Wo  - 
- (’o  ♦ f ■'»>  <•’1  * hd ) 


■o' 


Ik 


« ■ STTr  Ui  ‘i  1 S;  [‘i  J ■’i 

11  o * 


T 


V = - — ^— - {[y  + - K 1 [Ji  + ij  - [Yi  + - kJ  [J  - I ]} 

J,  + I,  ‘-  O TT  0-'  ^ 1 1-*  ^1  IT  P *-  O o-' 


'^l-l-'l  "l' 


where  the  arguments  of  the  Bessel's  functions  are  aa  = (— 


The  parameters,  R,  y and  o)^  are  independent  of  w and  represent 
a reflection  coefficient,  a mass  ratio  and  a characteristic  frequency  of  the 
plate,  respectively.  They  are  defined  by 


pet 


pet 

s s s 


^ 8y  a^ 
P 


li).  = (D  /y„)-"'Va 


1/2/2 


It  is  to  be  noted  that  is  not  the  fundamental  frequency  of 
the  plate  but  that  the  natural  frequencies  of  the  plate  itself  can  be  expressed 
in  terms  of  o)..  The  quantity  c is  a fictitious  plate  speed  defined  by  c = aoj 


3.3.  RESULTS  AND  NUMERICAL  EXAMPLES 

It  is  convenient  at  this  point  to  give  the  results  for  the  three 
input  cases  separately.  The  amplification  factor  for  the  case  of  ground 
shock  is  taken  as  ju/A^I  since  this  represents  both  the  ratio  of  displace- 
ments and  the  ratio  of  the  equipment  acceleration  to  the  ground  shock  accelera- 
tion. In  the  other  two  cases  the  amplification  factor  should  Involve  tlie 
equipment  acceleration  and  the  force  terms  and  these  are  of  different  dimen-- 
slon.  The  results  are: 

a.  Ground  Shock  Input 


s|  [N^  + R=  (~  ) Q^]- 


(3-27) 


Air  Shock  Input 


U/t 

= 


{[SN  + R^  (-  ) T(i]^  + R^  (-  ) [SQ  - TN]^}^ 


(3-28) 


where 


where 


p = — ^ — 2 is  a dimensionless  pressure. 
^ U„t  0)* 

P P * 


c . Shaker  Input 


U/a 


{[VN  + R^  (-  ) Zq]^  + R^  (-  ) [VQ  - ZN-J^' 


+ R^  (-  ) 


is  a dimensionless  force. 


(3-29) 


The  parameter  R expresses  the  reflection-transmission  characteristics 
of  the  interface  between  the  plate  and  the  sliell.  When  a pulse  travels  up  the 
shell  to  the  interface  some  energy  is  reflected  back  into  the  shell  and  some  is 
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transmitted  into  the  plate.  If  R is  less  than  1 the  reflected  wave  is  an 
unloading  pulse  and  R greater  than  1 means  that  a loading  wave  is  reflected 
back  into  the  shell.  The  coefficient  depends  on  the  relative  dynamic  stiff- 
ness of  the  two  elements  and  in  this  case  it  seems  likely  that  the  plate  be- 
ing a flexural  element  will  be  much  less  stiff  than  the  shell. 

The  existence  of  radiation  damping  in  the  system  is  due  to 
finite  values  of  R;  as  R ->  0 the  shell  becomes  infinitely  stiff  and  the 
equations  reduce  to  those  of  a clamped  plate,  and  as  R ->■  <»  to  those  of  a 
free  plate  (with  edges  prevented  from  rotating).  When  R ->■  0 each  of  the 
three  solutions  reduces  to  terms  having  only  N as  the  denominator.  The 
amplification  factors  thus  tend  to  infinity  at  the  zeroes  of  N.  In  the 
case  Y = 0 the  zeroes  of  N are  given  by 

(/uTw^)  I^  (/w/w*)  + (/(jj/w*)  (/oj/w^)  = 0 (3-30) 

which  is  the  characteristic  equation  for  the  natural  frequencies  of  a clamped 
plate.  The  roots  of  this  equation  are  tabulated  in  Jahnke  and  Emde  (Refer- 
ence 13).  When  R each  of  the  three  solutions  reduce  to  expressions  with 
Q as  the  denominator  and  for  y 0 the  zeroes  of  Q are  given  by 

(/u/w^)  I^  (i/w/u)^)  = 0 

which  is  recognized  as  the  equation  for  the  natural  frequencies  of  a free 
plate  with  edges  constrained  against  rotation. 

Based  on  the  dimensions  of  a typical  silo,  cases  have  been  con- 
sidered where  R = 0.01,  0.05  and  Y = 0.01,  .001,  .001  and  Equation  (3-27)  and 
have  been  computed.  The  results  for  the  air  shock  turn  out  to  be  very  similar 
to  the  ground  shock  case  and  have  not  been  included.  In  these  examples  it 
has  been  assumed  that  the  equipment  is  tuned  to  the  fourth  natural  frequency 
of  the  plate.  Typical  results  are  shown  in  Figures  8 and  9. 

It  is  clear  that  the  results  for  both  cases  are  qualitatively 
similar  to  the  earlier  results  on  the  more  simplified  model.  To  gain  insight 
into  the  origin  of  the  double  peak  for  this  more  complex  system  some  of  the 
asymptotic  properties  of  the  Bessel  functions  may  be  utilized  since  the  value 
of  (w/co^)  near  the  double  root  is  large. 


For  small  R the  amplification  peaks  are  given  by  the  zeroes  of 
N,  namely  by  those  values  of  cj  satisfying 
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(3-31) 


For  large  values  of  the  argument  of  the  Bessel  functions,  namely  /co/w^  (here 
Wq  = (12.58)^  co^,  12.58  being  the  fourth  root  of  the  eigenvalue  equation),  the 
Bessel  functions  can  be  approximated  by  asymptotic  expansions.  These  are: 

■^o  - J)  + lx 

^o  - f>  - lx 

-^1  = Vix  (X  - 1^  ) - 1^  sin  (X  - f^)} 

= VS  ^ 
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I (x)  = 
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I (s)  = 

Ki  (X)  = 

After  substituting  these  asymptotic  forms  into  Equation  (3-31)  and 
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This  expression  may  be  expected  to  be  highly  accurate  for  all  values  of 
the  argument  larger  than  the  lowest  frequency  of  the  plate  and  even  in  that 
case  the  accuracy  is  quite  good,  e.g.,  the  first  root  of  the  exact  expres- 
sion (Equation  3-30)  is  3.2  while  the  asymptotic  expansion  of  it  gives  3.14. 

The  roots  of  the  above  equation  are  given  by 
tan  (/w/to^) 


w/w. 


-Y 
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Figure  10  shows 


tan  (/oj/o)^^) 
to/w* 


-Y 


and 


— , for  the  case  where  the  equipment  frequency  to  is  tuned  to  a 

1 - to  /to  o 


natural  frequency  of  the  plate.  When  this  is  a high  mode  the  term  ('o/to*)  in 
the  denominator  is  very  large  and  tends  to  flatten  the  tangent  curve.  For 
fixed  nonzero  y the  roots  in  the  high  frequency  case  tend  to  be  located  at 
values  of  (to/to,,^  which  are  + •rT/2  distance  from  the  particular  natural 
frequency  of  the  plate  to  which  the  equipment  is  tuned.  For  low  frequency 
cases,  however,  the  roots  tend  to  coalesce  as  Y becomes  small.  There  are 
two  approximations.  In  the  low  frequency  case  the  spread  between  the  peaks 
of  the  double  resonance  is  given  by 

^ = 2 /T-y 
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and  for  high  frequency  cases 
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where  n is  the  mode  number  of  the  plate  natural  frequency  and  is  the  nth 

root  of  Equation  (3-30)  which  is  related  to  the  frequency  of  the  nth  plate 


mode  u)  through 
n 


w = A 
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These  results  indicate  the  spread  between  the  double  peaks  increases  as  the 
equipment  is  tuned  to  higher  plate  frequencies  for  low  order  modes  and  sub- 
sequently decreases  for  higher  order  modes.  The  phenomenon  of  the  double 
peak  is  thus  an  intermediate  frequency  one,  possibly  not  as  important  in 
either  the  low  range  or  the  very  high  range,  since  the  existence  of  the 
double  peak  provides  a range  of  frequency  over  which  the  system  is  suscep- 
tible to  resonant  conditions. 

Numerical  computations  of  the  amplification  factor  given  in 
Equation  (3-29)  are  shown  in  the  sequence  Figures  11a  through  lie  in  which 
results  for  shaker  excitation  are  shown  when  the  equipment  is  tuned  to  the 
first  five  natural  frequencies  of  the  plate  in  turn.  In  these  diagrams  the 
frequency  is  scaled  in  such  a way  that  the  equipment  frequency  is  unity  so 
that  tliere  is  a compression  of  the  frequency  scale;  in  terms  of  real  fre- 
quency the  increased  spread  in  the  double  peak  is  clear. 

Figure  12  shows  how  the  spread  varies  with  mode  number  of  the 
support.  The  spread  is  greatest,  offering  the  greatest  target  for  resonance, 
at  an  intermediate  frequency;  for  the  parameters  of  the  system  shown  in 
Figure  12,  tlie  fourth  mode  offers  the  largest  target.  Figure  13  shows  how 
the  amplification  of  the  equipment  at  resonance  varies  with  mode  number  of 
the  support.  Both  Figures  12  and  13  indicate  that  the  double  resonance 
phenomenon  is  critical  at  intermediate  frequencies.  Figure  13  shows  that, 
for  the  parameters  considered,  the  equipment  acceleration  can  be  as  much  as 
5,000  times  greater  than  the  input  acceleration.  The  results  in  these 
figures  include  only  the  effects  of  radiation-type  damping  (energy  radiated 
into  the  supporting  structure)  and  will  be  reduced  (probably  slgnific  ntly) 
when  structural  damping  is  Introduced. 
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SECTION  4 


ALTERNATIVE  MODELS  OF  DAMPING  FOR  HIGH  FREQUENCY  ANALYSIS 
4.1.  MAXWELL  MODEL  OF  DAMPING 

An  alternative  simple  modeling  of  structural  damping  which  has 
the  same  degree  of  simplicity  as  the  usual  damping  model,  but  which  it  is 
believed  will  lead  to  an  improved  estimate  of  structural  response  at  high 
frequencies  is  introduced.  The  model  is  based  on  the  well  known  Maxwell 
model  of  viscoelasticity  and  thus,  in  a one  degree-of-f reedom  system,  re- 
quires only  one  damping  parameter. 

It  is  known  from  analyses  of  wave  propagation  in  viscoelastic 
materials  that  given  the  appropriate  type  of  viscoelastic  model  a continuum 
can  propagate  discontinuities  as  waves.  The  essential  requirement  is  some 
form  of  instantaneous  elasticity  in  the  model.  The  material  will  also  be 
dissipative  and  the  discontinuity  may  be  attenuated  as  it  propagates  through 
the  material,  but  the  presence  of  dissipation  by  itself  is  not  a counter 
to  the  existence  of  a discontinuity.  There  are  models  of  material  behavior, 
for  example  the  Navier-Stokes  theory  for  fluids  and  the  Kelvin-Voigt  theory 
for  solids,  which  do  not  permit  the  existence  of  discontinuities.  In  both 
these  models  there  can  be  no  acceleration  waves  and  essentially  this  result 
arises  from  the  lack  of  an  Instantaneous  elastic  response  in  these  models. 

The  notions  of  Instantaneous  response  and  acceleration  waves  in 
extended  continue  are  very  closely  related  to  high  frequency  behavior  in 
bounded  systems,  and  thus  it  is  obvious  that  much  the  same  considerations 
will  pertain  to  these  situations. 

In  structural  dynamics  damping  is  generally  Included  in  the 
structural  model  in  the  form  of  viscous  damping  or  the  modification  of  it 
which  is  referred  to  as  structural  or  hysteretic  damping.  The  viscous 
damping  structural  model  is  an  analogue  of  the  Kelvin-Voight  viscoelastic 
continuum  model,  and  thus  it  seems  clear  that  high  frequency  response  and 
short  time  response  would  be  better  modeled  by  the  use  of  a structural  model 
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which  is  analogous  to  the  Maxwell  viscoelastic  continuum  model.  Such  a 
model  is  developed  in  the  following  section  and  it  is  shown  that  it  needs 
no  more  information  to  specify  it  than  does  the  usual  viscous  damping  model 
and  that  many  of  the  characteristics  of  viscous  damping  carry  over  to  the 
new  model. 


Although  it  is  relatively  straight  forward  to  specify  the  t lastic 
and  inertial  characteristics  of  a structural  system  it  is  very  much  more 
difficult  to  specify  the  damping  characteristics.  Normally  a rule  of  thumb 
approximation  is  used  in  which  the  damping  in  each  model  is  assigned.  Exactly 
the  same  form  can  be  used  in  this  approach. 


4.2.  GOVERNING  EQUATIONS  OF  MAXWELL  MODEL 


The  most  effective  way  to  illustrate  the  behavior  of  the  damping 
model  suggested  here  is  to  study  the  response  of  a single  degree-of-f reedom 
system.  In  this  case,  the  equation  of  motion  is  governed  by 

mil  + F = P(t)  (4-1) 

with 
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To  express  these  equations  in  non-dimensional  variables,  define 
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where  p(x)  = P(t)/F^,  and  primes  denote  differentiation  with  respect  to  x. 
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For  the  first  example,  steady-state  response  under  sinusoidal  load- 

P F 

ing,  i.e.,  P(t)  = P^  sinwt  or  p(x)  = R sinBx,  R = o/  o,  is  B = investigated. 

By  eliminating  f from  Equation  (4.3) 

y'"  + 2^y"  + y'  = R (2C  sinSx  + B cosBx) 

and  the  steady-state  amplitude  y is  then  given  by 

V = R i /F(1-B^)^  TT2gB)^  (T+2^2B‘*)  + (2CB)^ 

^ 6^  (1-B^)^  + (2^B)^ 

for  5=0,  the  amplitude  for  an  undamped  system  is  recovered, 

R 

y (1-B^) 

Moreover,  for  the  undamped  case, 

lira  9 = R 
B->0 

i.e.,  as  the  frequency  of  the  load  approaches  zero,  the  displacement  is 
bounded  and  approaches 

P 

o 

“o  K 

o 

F 

K = 
o u 

o 

For  5 different  from  zero,  however  y grows  without  bound  as  B"*'0.  This  is 
indeed  expected,  since  the  Maxwell  model  behaves  like  a fluid  under  small 
displacement^rates . Figure  14  shows  a plot  of  y/R  versus  B where  this  is 
clearly  exhibited. 

4.3.  GENERALIZATION  TO  MULTI  DEGREE-OF- FREEDOM  SYSTEMS 

In  view  of  the  fact  that  the  most  common  method  of  introducing 
damping  In  structural  dynamics  problems  is  to  resort  to  the  Kelvin  model  of 
linear  viscoelasticity,  a comparison  of  Kelvin  and  Maxwell  models  is  in  order. 
Consider  the  dynamic  equations  of  equilibrium  for  a multi  degree-of-f reedom 
system; 


. ^ 


AO 


with 


M ii  + F = P 


F = K u - C*  F 


These  two  equations  lead  to  the  third  order  linear  differential  equation 

M ii  + C*  M ii  + K u = C*P  + P (4-5) 

hot  C = C*  M and  integrate  Equation  (4-5)  to  obtain 


Mii+Cu  + Ku  = 


P + C M ^ P dt  + (C  u + K u - f1 

On  the  other  hand,  assuming  that  the  Internal  forces  are  governed  by  a Kelvin 
mode  1 , i . e . , 

F = K u + C u 

the  displacement  equations  of  equilibrium  take  the  form: 

Mu+Cu+Ku=P  (4-6) 

Thus,  both  the  Kelvin  and  Maxwell  models  lead  to  the  same  differentia]  equation, 
but  have  different  forcing  functions.  This  implies  that  a computer  program 
designed  for  solving  Kelvin-type  problems  can  be  easily  modified  to  solve 
systems  with  Maxwell-type  damping.  In  fact,  this  can  be  accomplished  even 
in  the  pre-processing  stage  of  computing  the  loads  for  a system. 

A significant  difference  exists  between  Equation  (4-5)  and  Equa- 
tion (4-6)  in  initial  conditions  that  need  to  be  specified.  Equation  (4-6) 

requires  that  the  displacement  u and  velocity  u at  time  t be  known  to 
' ~o  ~o  o 

obtain  the  solution;  these  quantities  completely  determine  the  force  F at 

t via 
o 

F = K u + C u 
^ ~ ~o  ~ -o 

With  a Maxwell  model,  however,  these  quantities  do  not  determine  the  initial 

value  of  F;  in  fact,  the  value  of  F at  t must  be  specified  in  addition  to 
~ o 
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and  u^.  This  aspect  of  the  Maxwell  model  suggests  that  it  might  indeed 
be  more  suitable  than  the  Kelvin  model  for  Impact  problems.  For,  while  a 
structural  system  at  rest  acquires  an  initial  velocity  following  an  impact,  it 
is  not  reasonable  to  assume  that  this  also  results  in  non-zero  values  for  re- 
sisting forces.  The  use  of  a Kelvin  model  does  incorporate  such  an  assumption, 
and,  if  the  damping  matrix  C is  proportional  to  the  stiffness,  K it  leads  to 
high  accelerations  and  low  displacements  in  the  higher  modes.  In  contrast, 
the  Maxwell  model  results  in  more  flexible  response  and  gives  lower  accelera- 
tions and  higher  displacements.  With  mass  proportional  damping,  modal  damp- 
ing is  inversely  proportional  to  modal  frequency  and  both  models  give  essen- 
tially the  same  answers  for  the  low  damping  ratios  resulting  from  such  an  as- 
sumption. Of  course,  many  of  the  difficulties  attendant  on  the  use  of  viscous 
damping  in  structural  analysis  of  impact  problems  or  high  frequency  analysis 
could  be  alleviated  by  the  use  of  general  viscoelastic  models  to  represent 
the  response  of  the  system.  However  the  degree  of  complication  which  results 
would  be  considerable  for  a large  system  and  it  would  be  almost  impractical 
to  obtain  the  necessary  experimental  data  to  define  the  response.  What  is 
needed  is  some  model  which  has  the  simplicity  of  the  viscous  model  which 
can  be  •enerally  defined  simply  by  specifying  a damping  factor  for  each  mode. 

It  seems  likely  that  the  Maxwell  model  when  converted  to  the 
structural  problem  might  be  useful  for  high  frequency  analysis.  It  should  be 
noted  that  the  same  number  of  parameters  are  needed  to  define  the  Maxwell 
model  as  are  needed  for  the  Kelvin  model  and  that  the  use  of  modal  damping 
factors  would  also  suffice. 
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SECTION  5 


APPLICATION  OF  RESULTS  TO  E.D.C.  TEST 

This  theory  can  be  applied  to  the  forthcoming  series  of  tests  on 
the  Electrical  Distribution  Center  (EDC)  to  be  conducted  for  DNA  during  the 
Summer  of  1977.  The  equipment,  part  of  the  safeguard  system,  in  this  system 
is  carried  on  a rectangular  plate  suspended  by  four  shock  isolation  devices. 
Thus  the  system  is  close  to  the  theoretical  model  considered  in  Section  3 of 
this  report  and  it  is  possible  that  some  items  of  equipment  carried  on  the 
isolated  plate  may  in  fact  be  tuned  to  some  natural  frequency  of  the  system. 
Shaker  excitation  applied  to  the  roof  of  the  building  will  correspond  to  the 
ground  shock  solution  and  shaker  excitation  applied  to  the  plate  will  corres- 
pond to  the  shaker  excitation  problem.  Acceleration  measurements  on  various 
items  of  equipment  on  the  isolation  system  under  these  two  types  of  excitation 
would  be  useful  and  revealing. 

The  main  difficulty  in  making  quantitative  predictions  of  the 
response  of  the  EDC  during  the  testing  program  is  of  course  the  lack  of  de- 
tailed information  on  the  building,  on  the  platform  and  on  the  Isolators 
themselves.  Some  information  has  been  provided  by  Waterways  Experiment  Sta- 
tion (WES),  Vicksburg,  who  is  conducting  the  tests  on  the  building  and  plat- 
form which  has  enabled  estimates  to  be  made  of  the  platform  loading  and  its 
bending  stiffness.  On  the  basis  of  this  data  the  platform  has  been  modeled 
as  an  Isotropic  plate  with  an  effective  moment  of  inertia  of  100  in."*  per 
inch.  The  weight  of  the  electrical  equipment  carried  by  the  platform,  report- 
ed by  WES  as  60,000  pounds,  was  assumed  to  be  uniformly  distributed  result- 
ing in  a plate  mass  per  unit  area  (including  the  self  weight)  of  0.75  pounds 
per  square  inch. 

Information  from  Barry  Burbank,  the  manufacturers  of  the  isolators 
(References  14  and  15)  suggests  that  an  isolator  can  be  considered  as  a spring 
which  when  carrying  a mass  of  20,000  pounds  gives  a frequency  of  0.73H.  This 
leads  to  a spring  stiffness  k given  by 
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I of  1,100  pounds/inch  for  each  isolator. 

t 

j The  platform  plan  dimensions  are  very  nearly  square  with  average 

side  length  25  feet.  The  fundamental  mode  of  a completely  free  square  plate 
is  one  in  which  there  are  nodal  lines  through  the  center  of  the  plate  parallel 
to  the  sides.  The  exact  solution  for  this  case  is  not  known  but  an  approxlma- 
tion  to  the  frequency  can  be  obtained  using  Rayleigh's  principle  assuming  a 
mode  shape  of  the  form  w = xy.  In  the  particular  problem  considered  herein 
the  same  mode  shape  is  maintained  and  the  contribution  of  the  Isolators  to 
the  potential  energy  is  included.  Assuming  a steady  harmonic  motion  of  the 
form  w = 4w^/a^  xy  cos  U)t , the  strain  energy  stored  in  the  plate  is 

32E  h^  w ’ 

O 9 

V = — - cos^  tot 

P 3(1+V)  a^ 

where  E is  Young's  modulus,  h is  the  half  thickness,  a is  the  length  of  the 
edge,  V is  Poisson's  ratio  and  o)  is  the  circular  frequency  of  oscillation. 

For  the  parameter  values  utilized  here  this  yields  '''p  “ pounds  per 

inch,  where  E has  been  taken  to  be  30  x 10®  psl  and  v = 1/3.  The  total 
potential  energy  stored  in  the  four  isolators  is  given  by 

V,  = 4 • k w ^ cos^  tot 
I 2 o 

I which  results  in  = 2,200  w^^  pounds  per  inch.  Thus  the  total  potential 

f 2 

■ energy  of  this  subsystem  is  402,000  w pounds  per  incli.  The  kinetic  energy 

\ 

' of  the  plate  is  1 

T = p h a'*  w ^ 0)^  sin^  tot  I 

9 o 1 

where  p is  the  mass  per  unit  of  volume.  Applying  Rayleigh's  principle,  the 

fundamental  frequency  is  estimated  as  23  Hz.  ' 

The  next  highest  natural  frequency  is  one  in  which  the  diagonals 
of  the  plate  are  nodal  lines.  Rayleigh  (Reference  16)  estimates  that  this 


44 


frequency  will  be  about  30  percent  liigher  than  the  fundamental  one,  resulting 
in  a figure  of  30  Hz. 

These  calculations  indicate  that  the  double  peak  amplification 
phenomenon  could  show  up  in  tiie  response  of  tlie  platform  in  the  frequency  range 
20  Hz  to  35  Hz.  The  frequency  spread  at  a double  peak  dictates  the  sliifting 
down  of  the  frequency  range  lower  limit  and  shifting  up  of  tlie  upper  limit. 

It  is  difficult  to  estimate  the  spread  since  the  effective  mass  of  the  en- 
tire structure  is  not  known.  The  actual  mass  of  the  entire  structure  exceeds 
1,200,000  pounds  but  the  contribution  from  backfill  is  not  known.  The  total 
mass  of  the  platform  is  around  70,000  pounds  indicating  a mass  ratio  around 
0.00.  For  such  a mass  ratio  tlie  spread  is  about  20  percent  of  the  central 
frequency. 

From  information  provided  by  WES  the  structure  Itself  appears  to 
have  a number  of  resonances  below,  in  and  above  this  range  and  thus  if  one  of 
the  structural  frequencies  is  close  to  one  of  the  platform  frequencies  the 
double  peak  should  appear  in  the  response  of  the  platform. 

It  is  therefore  recommended  that  the  test  series  be  planned  in  such 
a way  that  this  phenomenon  can  be  observed  if  it  should  occur.  Acclerometers 
should  be  placed  on  the  platform  at  each  corner  and  at  the  center.  Accelero- 
meters should  also  be  located  at  various  points  on  the  building.  The  shaker 
should  be  located  on  the  building  roof.  The  results  of  excitation  in  the  fre- 
quency range  15  Hz  to  45  Hz  within  which  several  building  resonances  are  expect- 
ed will  be  of  primary  interest.  In  particular  in  tlie  range  20  Hz  to  35  Hz  the 
possibility  of  double  peak  amplification  due  to  tuning  is  anticipated.  If  tun- 
ing occurs  it  is  expected  tliat  the  spread  between  peaks  would  be  of  the  order 
of  5 Hz. 

Experimental  results  of  dynamic  tests  on  the  isolators  by  the  manu- 
facturers (Reference  17)  suggest  that  the  Isolators  are  subject  to  a subsystem 
resonance  at  around  50  Hz  and  tlius  might  be  expected  to  interact  with  a plat- 
form resonance  around  this  value  if  such  should  occur.  The  actual  magnitude 
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of  this  internal  resonance  will  depend  on  the  pendant  length  and  fixity  condi- 
tions and  may  vary  from  situation  to  situation. 

There  are  disadvantages  in  restricting  attention  to  these  measure- 
ments alone  since  the  vibrational  characteristics  of  the  shock-isolated  plat- 
form of  the  isolators  or  of  the  attached  equipment  are  not  accurately  known. 

To  reveal  the  existence  of  the  phenomenon  postulated  in  this  report  it  would 
be  necessary  to  add  a small  structure-equipment  system  to  the  platform  which 
would  undergo  motions  during  the  experiment.  The  properties  of  each  element 
of  this  system  would  be  known  much  more  accurately  than  those  of  the  existing 
utility  and  electronic  substations.  The  new  equipment  element  could  be  adjust- 
ed to  be  both  non-tuned  and  tuned  to  several  platform  frequencies.  For  this 
purpose,  a simple  and  inexpensive  mounting  would  be  adequate.  A flat  plate 
square  in  plan  supported  on  four  angles  at  the  corners  is  suggested.  On  this 
plate  to  simulate  equipment  would  be  a set  of  reed  gages  of  a variety  of 
tuned  and  non-tuned  systems.  Alternatively,  a small  number  of  specially  de- 
signed spring-mass  oscillators  with  accelerometers  may  be  used. 

The  results  sought  in  these  tests  would  be  the  determination  of 
whether  the  double  peak  resonance  occurs  and  if  so  at  the  frequencies  predicted. 
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SECTION  6 


FUTURE  RESEARCH  DIRECTIONS 

The  results  described  in  this  report  suggest  certain  directions 
for  future  research  in  the  analysis  of  equipment-structure  systems  at  high 
frequency.  This  work  falls  into  four  main  areas  with  substantial  interaction 
between  them. 

6.1.  DAMPING 

The  theory  developed  so  far  should  be  extended  to  include  struc- 
tural damping  which  is  expected  to  play  an  important  role  in  the  response  of 
the  system.  Previous  work  has  shown  that  conventional  structural  damping 
models  are  not  realistic  at  high  frequencies  and  a new  approach  to  damping 
suitable  for  high  frequencies  should  be  developed. 

6.2.  MULTI-MODE  RESPONSE  OF  ATTACHED  EQUIPMENT 

In  the  previous  work  the  equipment  element  of  the  structure- 
equipment  system  has  been  modeled  as  a single  degree-of-f reedom  system.  This 

should  be  expanded  to  models  which  include  multiple  degree-of-f reedom  behavior  ' j 

in  the  attachment.  This  may  reveal  important  coupling  phenomena  overlooked 
in  the  simple  model. 

6.3.  EXTENSION  TO  GENERAL  STRUCTURAL  SYSTEM 

It  has  become  clear  in  the  course  of  the  previous  work  on  the 
idealized  systems  that  many  of  the  results  are  not  special  to  these  models, 
but  will  appear  in  any  structural  system.  It  should  be  possible  to  develop 
these  results  for  more  general  structural  systems  utilizing  a technique  of 
cascading  Impedance  functions  for  subsystems.  These  Impedance  functions  may 
be  obtained  either  by  experiment  or  by  suitable  numerical  analysis.  Such 
information  could  serve  as  input  to  models  of  the  kind  used  here  to  provide 
the  response  of  real  structure-equipment  systems. 
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6.4.  TRANSIENT  RESPONSE 


Previous  work  has  been  concerned  with  steady  state  response. 
Techniques  needed  to  synthesize  these  results  for  application  to  transient 
loading  typical  of  air  shock  or  ground  shock  should  be  studied. 

The  scope  of  the  work  outlined  above  is  broad.  It  must  be  com- 
pleted before  work  can  begin  in  adapting  finite  element  techniques  to  high 
frequency  problems.  However,  a preliminary  study  of  integration  algorithms 
which  may  be  useful  in  treating  this  class  of  problems  should  be  initiated. 
If  the  work  described  above  goes  smoothly,  some  experimentation  with  an 
existing  finite  element  code  would  be  appropriate. 
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Figure  3.  Diagrammatic  Representation  of  Advanced  Model 


structure  impedance 
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Diagrammatic  Representation  of  Silo-Like  Structure  with  Attached  Equipment 
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Figure  9.  Amplification  Factor  as  a Function  of  Non-Dimensional  Frequency 
Shaker  Input 
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Figure  10.  Solution  Curves  for  Tuned  Resonance 


amplification  factor 


Figure  11a.  Amplification  Factor  as  a Function  of  Nondimensional  Frequency 
for  Shaker  Problem  with  Equipment  Tuned  to  Various  Plate 
Frequencies.  R = .01,  y = .001. 
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Figure  11 
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Amplification  Factor  as  a Function  of  Nondimensional  Frequencj 
for  Shaker  Problem  with  Equipment  Tuned  to  Various  Plate 
Frequencies.  R = .01,  y = .001. 


Equipment  Tuned  to  Third  Plate  Frequency 
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Figure  lid.  Amplification  Factor  as  a Function  of  Nondlmenslonal  Frequency 
for  Shaker  Problem  with  Equipment  Tuned  to  Various  Plate 
Frequencies.  R = .01,  y = .001. 


Equipment  Tuned  to  Fourth  Plate  Frequency. 


AMPL  I F I CAT  I ON 


2 5 


OMEGA 


Figure  He.  Amplification  Factor  as  a Function  of  Nondimenslonal  Frequency 
for  Shaker  Problem  with  F.qulpment  Tuned  to  Various  Plate 
Frequencies.  R = .01,  Y “ .001. 

F.qulpment  Tuned  to  Fifth  Plate  Frequency 
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MAXIMUM  AMPLIFICATION  FACTOR 


Figure  13.  Maximum  Amplitude  of  Double  Resonance  as  a Function  of  the 
Mode  Number  to  which  the  Equipment  is  Tuned 
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